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Population transfer between two identical, communicating defects in a one-dimensional tight- 
binding lattice can be systematically controlled by external time-periodic forcing. Employing a 
force with slowly changing amplitude, the time it takes to transfer a particle from one defect to 
the other can be altered over several orders of magnitude. An analytical expression is derived 
which shows how the forcing effectively changes the energy splitting between the defect states, and 
numerical model calculations illustrate the possibility of coherent control of the transfer. 

PACS numbers: 73.43.Jn, 71.23.An, 72.10.Fk 
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I. INTRODUCTION 

The possibility to control the state of a quantum sys- 
tem by external forcing, or even to create desired target 
states by judiciously designed electromagnetic pulses, is 
of substantial conceptual and practical interest. In this 
paper, a specific scenario is investigated which emerges 
in one-dimensional tight-binding lattices under the in- 
fluence of periodic forcing, such as laser-irradiated poly- 
mers, or semiconductor superlattices interacting with far- 
infrared radiation. If the on-site energy is deliberately 
modified in an identical manner at two sites of the other- 
wise regular lattice, two defect states emerge which "com- 
municate" if they are not too far apart, i.e., they possess 
localized eigenstates connecting both defects. It will be 
shown that the degree of communication, expressed by 
the splitting of the energies of the two defect states, can 
be varied over several orders of magnitude if the system 
is subjected to strong periodic forcing. From this fol- 
lows the possibility to achieve well-controlled population 
transfer from one defect to the other, if the forcing am- 
plitude is suitably shaped. 

The paper relies on the framework provided by quan- 
tum mechanical Floquet theory for periodically forced 
systems, which has repeatedly been found useful in 
the analysis of solid-state devices driven by external 
f orcea i i 2 J 3 J 4 The material is organized as follows: Secs.HTI 
and IIIII briefly provide the required background on local- 
ization at isolated defects, and on the influence of a peri- 
odic external force on a defect state. In Sec. llVl an analyt- 
ical expression for the energy splitting between two defect 
states is derived, and generalized to describe the relevant 
quasienergy splitting when the periodic force is turned 
on. Sec. [V] explains in detail the strategy for achieving 
controlled population transfer between the two defects. 
Conclusions are drawn in the final Sec. I VII 



II. LOCALIZATION AT A SINGLE DEFECT 

A single particle on a one-dimensional, infinite, regu- 
lar lattice with matrix elements connecting neighboring 
sites only, as is appropriate in the tight-binding limit, is 



described by the Hamiltonian 

H = 



— E {|*x*+ii + i*+i><*|}i w 



where \£) is the Wannier state localized at the l-th site, 
adopting the normalization (£'{£) = 5j>i,£. Denoting the 
lattice constant by d, its eigenstates 
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are extended Bloch waves, labeled by the wave number k. 
The hopping matrix elements — W/4 in Eq. have been 
chosen such that the energy dispersion reads 



W 

E(k) = -— cos(kd). 



(3) 



corresponding to a band of width \ W\. For positive W , its 
minimum lies at k = 0. In order to introduce a defect into 
the ideal system Q , the on-site energy at the site 7 is now 
altered by an amount v, giving rise to the perturbation 



^■ = |7M7l, 7>0. 
The single-defect Hamiltonian 

H = H + V t 
then admits a localized state with energy 

TF 



E =V 



where 



4^ 2 
W 1 



— > 



w 



< 



(4) 
(5) 

(6) 

(7) 



so that E falls either above (for v/W > 0) or below 
(for v/W < 0) the energy band The probability to 
find the particle at the ^-th site, when it is bound by the 
defect, is given by 
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as can be derived with the help of resolvent operator tech- 
niques (see, e.g., Ref. 5). In the following, knowledge of 
the amplitudes ae themselves will be required: As shown 
in the Appendix 1X1 the eigenstate corresponding to the 
defect energy |(BJ| reads 



with 
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where Af is the normalization constant, and the auxiliary 
quantities 



x± 



Av 2 

w 2 



1±2 
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w 



(ii) 



have been introduced. They obey x-x+ = 1, and 
x- < X+. The fact that the wave function JJJJ with 
the amplitudes IjlUfl indeed is normalizable for ^ ^ 0, 
and hence describes a localized state, follows immediately 
from < |x_ I < 1. 



III. LOCALIZATION CONTROL BY FORCING 

If an oscillating electric field is applied to the system, 
linearly polarized along the direction of the lattice and 
with amplitude F, the interaction is modeled by 



H int (t) = eFdcos(ujt) \() t {(] 



(12) 



where e is the particle's charge. Since the total Hamil- 
tonian H + V T + H- m t(t) then is periodic in time, with 
period T = 2tt/u>, the Floquet theoremSiSi asserts that 
there is a complete system of wave functions of the form 



driven particleii, an effect which should be observable in 
far-infrared driven semiconductor superlatticesi^. 

Without periodic forcing, the strength of a defect of 
the form is determined not by the on-site energy v 
alone, but rather by the ratio v/W: The larger |i//W|, 
the shorter is the localization length of the defect state, as 
witnessed by Eq. JHJ. This leads to the conjecture that in 
the presence of periodic forcing the defect Floquet state 
again is described by Eq. JSJ, but with v/W replaced by 
v/Weft, so that the localization length of the defect state 
should become strongly dependent on the amplitude of 
the forcing. In particular, when eFd/(huj) approaches a 
zero of Jo, the state should be confined entirely to the 
defect site, if only v £ 0. 

As shown in Ref. |5j, this conjecture indeed is correct 
in the high-frequency regime, where Hlj is significantly 
larger than the bare band width \ W\. This is illustrated 
in Fig. for a defect with strength v/W = 0.1: The up- 
per panel shows the occupation probabilities p n = \a n \ 2 
for a defect located at the site n = in the absence 
of the periodic force; the discrete values have been con- 
nected by lines to guide the eye. The lower panel 
shows the state under the influence of a periodic force 
with high frequency hu>/W = 7.5 and scaled amplitude 
eFd/(huj) = jo,i — 2.4048, equal to the first zero of J : 
As expected, the state now is confined almost entirely to 
the defect site. Hence, the extension of the defect state 
is governed by the driving force's amplitude. This ef- 
fect will be exploited in Sec.^to control the population 
transfer between two communicating defects. 



IV. ENERGY SPLITTING AND QUASIENERGY 
SPLITTING FOR COMMUNICATING DEFECTS 



In the following, the real amplitudes ai for the state 
bound by a single defect placed at a site 7 > are nor- 
malized such that 



m))=e- iet/H \u(t)) , 



(13) 



where \u(t)} is a T-periodic function. The quantity e is 
called "quasienergy" , in analogy to the quasimomcntum 
in solid state physics. 

If there is no defect, the quasienergies for a periodically 
driven particle in a tight-binding lattice with nearest- 
neighbor coupling readSiifi 

W ( eFd\ 

e(k) = J costkd) mod huj , (14) 

2 \ Tvjj ) 

so that the bandwidth is effectively quenched according 
to the zero-order Bessel function Jq: 



n-.,r = n-. M '^) 



(15) 



At the zeros of Jo, the band "collapses". This band 
collapse manifests itself as dynamical localization of the 



(16) 



which implies that the normalization constant M in 
Eq. I|10|) is now given by 



7 



-1/2 



E 



2(-y-») + x 2 ^ 7 ' 1 

i—1 i="y-j-l 



1 - x 2 



l + x 2 _- x 2 J 



(17) 



This defect state obeys the Schrodinger equation 

'60 + Vr) Ho) - E |Vo> • (18) 



Next, a second, identical defect is introduced into the left 
half of the lattice at the site labeled —7, as described by 
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FIG. 1: Occupation probabilities p n for a defect state in a 
one-dimensional lattice Q with a defect of the form @ at 
the site n — 0; the defect strength is v/W = 0.1. The upper 
figure shows the state in the absence of periodic forcing; it is 
exponentially localized according to Eq. JSJ. The lower figure 
shows the defect Floquet state when the system is driven with 
frequency hw/W = 7.5 and amplitude eFd/(huj) = 2.4048, 
corresponding to the first zero of Jo. Since the effective band 
width 1151 is zero here, the defect state is localized almost 
entirely at the defect site. 



V\ = |-7>"(-7l i 7>0 



(19) 



The two defects then carry two localized states |V>i) and 
|^2), obeying the eigenvalue equations 



Ho + V T + Vl) |Vi) = E x |Vi) , (20) 
H + Vr + Vi) IV2) = E 2 |V> 2 ) • (21) 



By symmetry, good approximations to these two defect 
states are given by the even and odd linear combinations 



KM = E 7fi( ai± a - l) M ■ 

£— — 00 

From these, auxiliary wave functions 
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x - 
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1=1 v 2 



(22) 



(23) 
(24) 



are defined, which have nonvanishing amplitudes in the 
right half of the lattice only. Forming the scalar product 

of Eq. JTHJ with (ipi | then gives 

fa\H a + Vr |Vo) - E ^h\ ifo) 5 (25) 
forming that of Eq. (|20l) with (ip 1 leads to 

<^\H + V r = V>0 , (26) 

since (ipo\Vi = (see Eqs. ((TUJ and (ESI)- By definition, 

(V'll V'o) = ("00 1 V>i)- Hence, subtracting the above two 
equations yields 

(E Q - £1) Wi\ lfo> = W1W0 |Vo) - Wv\Ho |V>1> ■ (27) 

It is now stipulated that the localization of the state \tpo) 
around the site 7 > be sufficiently strong that its am- 
plitudes in the left half of the lattice are negligible. One 
then has 



(V'll tpo) ^ -^(^ol $0) = 77= 



(28) 



leading to 

Eq-E\ ~ --^ | (ai+a_i) a e+1 + (ai+q_i) ai-ij 
VF 



00 



+ — X! ( a <?+i+ a -«-i) 
1=1 

+d£ (a £ _i+a_ £+ i) I 

-— ^ j(a_^+i - a_ t+1 at) 
1=1 

.)}■ 



+ a- 



a-e-iae) 



(29) 



These telescope series can be summed immediately 
(En=i _ = b N - b and lim^oo a £ = 0), re- 

sulting in 



Eq — Ei 



W 



(ai - a-i)a 



In the same manner, one also derives 

W 

E 2 ~ E = — (ax - a_i) a 



(30) 



(31) 



Summing these two equations finally leads to a surpris- 
ingly simple expression for the splitting of the energies 
associated with the two defects: 



Ei — E\ 



W 



(ai - a-i)a 



(32) 



Since (a\ — a_i) /2 can be taken as the discrete derivative 
of the wave function at the site £ = 0, the energy split- 
ting is determined by the product of the wave function 
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FIG. 2: Relative splitting of the quasienergies in a periodically 
driven lattice with two identical defects of strength u/W — 0.1 
situated at 7 = ±3, under the influence of forcing with fre- 
quency hui/W = 7.5. The full line shows the exact, numeri- 
cally obtained data. The dashed line results from the Herring- 
type approximation 1341 , with the bandwidth W replaced by 
the effective band width 1151 . 



itself and its derivative halfway between the two defects, 
i.e., by the current prevailing there. Thus, this Eq. I|32(l 
constitutes a discrete analog of Herring's formula, which 
describes the tunneling splitting for wave functions in 
double well potentialsi^i^i^. In view of Eqs. I|1(J|) and 
(|17|l . one then has 
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As with a single defect, the energy eigenvalues are re- 
placed by the corresponding quasienergies in the presence 
of time-periodic forcing, and the energy splitting turns 
into a quasienergy splitting. The above result , with 
W replaced by W e g according to Eq. l(T5|) , should also be 
a good approximation to the quasienergy splitting if the 
driving frequency is sufficiently high. 

In order to check this hypothesis, the time-dependent 
Schrodinger equation for the periodically forced two- 
defect system has been solved numerically, and the 
quasienergies for the localized states have been obtained. 
The results for two defects at 7 = ±3, with v/W = 0.1 
and hjjjW = 7.5, are plotted in Fig.[21as functions of the 
scaled amplitude eFdj (hw). As can be seen, the agree- 
ment between the analytical approximation and the exact 
numerical data becomes excellent when eFd/(huj) > 1. 



V. CONTROLLED POPULATION TRANSFER 

It is assumed now that initially, at time t = 0, the 
particle is localized at one of the two defects, and the 



periodic force is present. Denoting the two Floquet func- 
tions associated with the defects at the given driving am- 
plitude F by \uf^ 2 (t)), and their quasienergies by ef 2 , the 
initial state is given by a superposition 



1 



V2 



(|uf(0))±|uf(0)>) 



(35) 



Under the influence of forcing with constant amplitude, 
this state evolves in time according to 



W)) 



L( luf(i))±|uf , (t))axp(-iAe Ji 't//i)) 



V2V 
x exp(— ief t/h) 



where 



a F F F 

Ae* = e t 2 -e{ 



(36) 



(37) 



denotes the quasienergy splitting. Hence, the particle is 
coherently oscillating between the two defects; the trans- 
fer time T trans , after which the particle will be found at 
the other defect, is given by 



A^ 



(38) 



and thus depends on the driving amplitude F. 

When the amplitude changes sufficiently slowly in 
time, the system responds in an adiabatic manner—. 
Hence, under the influence of a slowly varying ampli- 
tude F(t), the initial state (|35|) follows the instantaneous 
Floquet states and evolves into 



W)) 
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V2 



f <*>(*)) 

uf t) (t))exp^i^ t drAe^) 



± 



x cxp j — ^ j dr ef ^ 



(39) 



This implies that the transfer time from one defect to the 
other now is given by the relation 



drAe F M 



(40) 



which constitutes an immediate analog of the 7r-pulsc- 
condition known from two-level systemsit This is the 
physics which will now be exploited for coherent control 
of population transfer between two communicating de- 
fects. 

To this end, the driving amplitude F(t) is shaped such 
that one has eFdj (fno) — j'0,1 for t < 0, so that the defect 
Floquet states are confined to their respective sites, and 
the communication between the two defects is effectively 
disrupted. Then the amplitude is adiabatically lowered 
such that the defect states start to overlap significantly, 
and the particle initially tied to one defect oscillates to 
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the other. If the amplitude then rises again and reaches 
the "collapse" value 

huj 

-Fcollapsc = jo,l"^ (41) 

at t = Tp U i sc , and is kept constant thereafter, the particle 
has been transferred to the final state, and will stay there. 
The time T pu i se is chosen such that 

- / drAe F(T ) =tp ; (42) 

to obtain a transfer from one defect to the other, 
again tp = it has to be chosen. 

For a matter-of-principle demonstration of this sce- 
nario, one may employ the envelope function 

F(t) = F co iiap Se (43) 

exp (- 2i£^) + a + cxp (- ^S^'H 

x 7 \ 

a + 1 + exp (-2^7 J 

for < t < Tpuiso, where T ram p quantifies the character- 
istic time interval during which the amplitude is ramped 
down and up again, with T -C T ramp <C T pu i se being un- 
derstood. The parameter a has been introduced in order 
to allow for a nonvanishing amplitude at intermediate 
times, at t « T pu i S e/2. 




Fig. 0] demonstrates that the desired result is obtained 
when the time scales are prolonged: With T ramp = 200 T 
and T pu i S c = 1070 T, one has practically complete trans- 
fer. 



°- 0. 




FIG. 4: Adiabatic population transfer at communicating de- 
fects: Now Tramp = 200 T, and T pu i ae = 1070 T; the other 
parameters are as in Fig. [3] 

It is, of course, also possible to employ the method 
outlined here to prepare the particle in superpositions 
of defect states. For instance, if one chooses T pu i se such 
that the phase integral in Eq. I|4U|) yields 7r/2, rather than 
7r, the resulting state describes a particle which, after 
initially being localized at a single defect, is eventually 
found with equal probability on cither one. In the same 
manner, any desired probability ratio can be obtained; 
Fig. [S] shows an example where the particle remains with 
a probability of 1/4 at the initial defect, and is found 
with a probability of 3/4 at the other one, after forcing 
according to Eq. 1(33} . 



VI. CONCLUSION 

It has been shown in this paper that a periodic force 
can drastically alter the energy splitting associated with 



FIG. 3: Incomplete population transfer resulting from fast, 
nonadiabatic changes of the amplitude function. There are 
two defects at 7 — ±3 with strength v/W = 0.1; the forcing 
frequency is huo/W = 7.5. The envelope function is given by 
Eq. with a = 1.5, T ramp = 50 T, and T pu i sc = 595 T. 

Plotted are the probabilities to find the particle at the two 
defect sites. 




The above considerations rely heavily on the adiabatic 
principle; if the parameter variation proceeds too fast, 
complete population transfer is not achieved. This is il- 
lustrated in Fig. O for a system with defect parameters 
7 = ±3 and v/W — 0.1, subjected to forcing with fre- 
quency hiv/W = 7.5 and the envelope function H43|l. set- 
ting a — 1.5, T r amp = 50 T, and T pu i sc = 595 T. Under 
these conditions, the transfer remains incomplete; slight 
oscillations visible in the occupation probabilities of the 
defect sites indicate non-adiabatic dynamics. However, 



1 1 

462.5 925 

t/T 

FIG. 5: Probabilities to find the particle at the two defects 
for T ram p = 200 T and T pu i sc = 925 T; the other parameters 
are as in Fig. This choice of parameters aims at producing 
a coherent, Schrodinger cat-like state, with probabilities of 
25% and 75% for finding the particle at the defects. Arbitrary 
other final probability ratios can be achieved as well. 
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the two states bound by two identical defects in a one- 
dimensional tight-binding lattice. In the presence of 
the force, the band width W entering the energy split- 
ting H34fl has to be replaced by the effective width (|15|l . so 
that the splitting can be monitored within wide ranges, 
and even be completely suppressed. Thus, the times re- 
quired for coherent population exchange between the de- 
fects can be varied over several orders of magnitude by 
adjusting the amplitude of the force. 

The strategy employed here to give a matter-of- 
principle illustration of coherent control of population 
transfer relies on the adiabatic principle, and thus is 
restricted to forces with slowly varying amplitudes. It 
appears possible to overcome this restriction: Utiliz- 
ing techniques developed for the coherent control of 
moleculesiSiiS, it might be possible to design even rapidly 
changing envelopes which effectuate a guided transport 
of a particle from one defect to the other, or to create su- 
perposition states with well-defined weights. This might 
open up new perspectives for the design of quantum log- 
ical devices. 
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Singling out the defect site, one has 
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pW\ — \a 1 - — (a 1 ^ l +a. l+l ) ) | 7 ) 
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Utilizing the explicit expression (|10|l for the ampli- 
tudes at, this leads to 
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APPENDIX A: EIGENFUNCTION FOR A 
SINGLE DEFECT 

In order to demonstrate that | -0o } as defined by 
Eqs. © and (fTTTfl indeed is an eigenfunction of the single- 
defect Hamiltonian © with the energy eigenvalue JHJl, 
it is helpful to invoke the definition Q, and to write 
Wp J jjp J . One then obtains 



\l)pW 



w\ 



w 



E iWit+M 



P w 



e=-c 
v 

W 



(7l^o) (Al) 

\t)(l-l\}\1k) 

a i+ i)\S>) . 



T E ^ 



According to the definition one has x_ l = x + . 

Hence, xZ 1 + — x + + a;_ = 2\J + 1, which im- 
plies 



W Uv 2 ^ 



(A4) 



This proves the assertion. 
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